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In this paper, we present the results of a study made for the effect of nuclear medium in the charged 
current induced quasielastic lepton production(CCQE) and the incoherent and coherent one pion 
production (CCl7r + )processes from 12 C in the energy region of 0.4-3GeV. The theoretical results 
are compared with the recent experimental results for the ratio of charged current induced one 
pion production cross section to the quasielastic lepton production cross section reported by K2K 
collaboration. We also present the results for the angular and momentum distributions of leptons 
and pions produced in these processes. 
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I. INTRODUCTION 



Recently K2K collaboration has reported experimental results on the inclusive single pion production induced by 
charged current by neutrinos in the energy region of 0.4-3 GeV[l[. Similar results were earlier reported by MiniBooNE 
collaboration in the energy region of lGeV 0. The study of the energy dependence of the pion production cross 
sections in neutrino reactions for nuclei in this energy is important in modeling the neutrino nucleus cross sections 
for various Monte Carlo neutrino event generators used in analyzing the present neutrino oscillation experiment at 
MiniBooNE, K2K and future experiments to be done by T2KQ and NO^A 0] collaborations. 

The present experimental result by K2K collaboration reports the energy dependence of the ratio of inclusive 
cross sections of the charge current induced pion production and quasielastic reaction in neutrino reaction on 
polystyrene ( C's H %) target. Both of these reactions are influenced by the nuclear medium effects in the energy re- 
gion of this experiment. The Monte Carlo simulation of single pion production in this reaction uses the Rein and 
Sehgal model [5| which does not include any nuclear medium effect in the production process. The NEUT generator 
used in the analysis of K2K experiment [lj], however, includes the nuclear effects arising due to final state interactions of 
pions with the nucleus like pion absorption and pion scattering |6] . The quasielastic reaction, on the other hand, uses 
Smith-Moniz model 0] which does not include the effect of nuclear medium arising due to nucleon-nucleon correlations 
but includes only the effect of Pauli principle and Fermi motion in a Fermi gas model. There are many calculations 
for the nuclear medium effects in the quasielastic reaction Q- 
pion production which take into account explicitly the nuclear medium and final state interaction effects in various 
energy region (26|- 38] . In this paper we present a study of nuclear medium effects in these reactions in the energy 
region of 0.4-3 GeV averaged over the neutrino flux of the K2K experiment. The single pion production in this region 
is dominated by the resonance production in which a A resonance is excited and decays subsequently to a pion and a 
nucleoli. When this process takes place inside the nucleus then there are two possibilities, the target nucleus remains 
in the ground state leading to coherent production of pions or is excited and/or broken up leading to incoherent pro- 
duction of pions. In the present paper, we have considered both the production processes in the A resonance model in 
the local density approximation to calculate single pion production from nuclei. The effect of nuclear medium on the 
production of A is treated by including the modification of A properties in the medium. Once pions are produced, 
they undergo final state interactions with the final nuclei. For the incoherent pion production final state interaction 
has been treated using a Monte Carlo code for pion nucleus interaction [39f . while for the coherent case, the distortion 
of the pion plane wave is calculated by using the Eikonal approximation [40] . In the case of quasielastic reaction, the 
effects of Pauli principle and Fermi motion are included through the Lindhard function calculated in local density 
approximation The details of these calculations are given in earlier publications [l9| , [26j , [37j • 

In this paper we report the results for the nuclear medium effects on total cross section <j{E), Q 2 distribution 
■jq? and momentum distribution ^ for the muon, in the case of quasielastic and single pion production induced 

a cclw+ (E) 

by charged currents. We compare our results of the ratio R(E) = ^ccqe jW ; with the recent observations of K2K 
collaboration. The results for the pion momentum and angular distributions have also been shown for the incoherent 
pion production process. In Sees. II & III, we briefly describe the various formula used in the present calculations for 
quasielastic and inelastic reactions and present the results and discussion in Sec. IV. 
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II. QUASIELASTIC REACTION 

The basic reaction for the quasielastic process is a neutrino interacting with a neutron inside the nucleus i.e. 

u^k)+n(p)^ ti-(k')+p(p') (1) 

The cross section for quasi-clastic charged lepton production is calculated in the local density approximation by 
taking into account the Fermi motion and the Pauli blocking effects through the imaginary part of the Lindhard 
function for the particle hole excitations in the nuclear medium. The renormalization of the weak transition strengths 
are calculated in the Random Phase Approximation(RPA) through the interaction of the p-h excitations as they 
propagate in the nuclear medium using a nucleon-nucleon potential described by pion and rho exchanges [4l| . The 
effect of the Coulomb distortion of muon in the field of final nucleus is also taken into account using a local version of 
the modified effective momentum approximation(MEMA) [19(. The total cross section a(E v ) for the charged current 
neutrino induced reaction on a nucleon inside the nucleus in a local Fermi gas model is written as 

9C v 2 ens 2 f) r r m^ rk'max /•! i 

a(E v ) = c / r 2 dr k' 2 dk' dcos6——L^J^ pA ImU N (q ,q), (2) 

71" Jr min Jk' min J-l 

where L^ v — ^ L^L^ and J^p^ = 2 S J^J U \ calculated with RPA correlations in nuclei. 
The leptonic current and the hadronic current are given by 

L^ = ii(k')^(l-j 5 )u(k) (3) 



J» = uip'^q 2 )^ + F 2 {q 2 )ia^^ + F A (q 2 )^ l5 + F P (q 2 )q^ 5 ]u(p)- (4) 

where q(— k — k') is the four momentum transfer, M is the mass of the nucleon, Gf(= 1.16637 x 10~ 5 GeV~ 2 ) is the 
Fermi coupling constant and 8 is the lepton angle. Un is the Lindhard function for the particle hole excitation [Til ] . 



7 
6 
5 

§4 
% 3 
2 
1 




; 1 1 1 1 1 1 1 1 


i i i i 1 i i i i 1 i i i i 1 i i 




~ r Incoherent 


— 


: WithW<2GeV 


E 


✓ 

— * / 


s ~. — ' 

s 

s ^> 




/ / 

: // 


■■ Free 




•' '/ 
— • // 


- With Medium effects(ME) 




- ■' y 

~- , *>r, , , , i 


- With ME + Pion Abspn 




0.5 1 


1.5 2 2.5 






E y (GeV) 





FIG. 3: Total scattering cross section for the induced incoherent charged current one pion production process in 1 C. 
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FIG. 4: Total scattering cross section for induced coherent charged current one pion production process in 1 C. 

The form factors Fi(Q 2 ), -F^Q 2 ), Fa(Q 2 ) and Fp(Q 2 ) are isovector electroweak form factors. For our numerical 
calculations we have used the parameterization of Bradford et al. [52] for the vector form factors Fi(Q 2 ) and F^Q 2 ) 
with dipole mass My=0.84GeV and a dipole form for Fa(Q 2 ) with dipole mass MA=l.lGeV. 

Inside the nucleus, the Q-value of the reaction and Coulomb distortion of outgoing lepton are taken into account 
by modifying the Lindhard function t/jv(9o,q) by [/jv(go — V c (r) — Q,l), where V c (r) is the Coulomb potential of 
the final nucleus. Furthermore, the renormalization of weak transition strength in the nuclear medium in a random 
phase approximation(RPA) is taken into account by considering the propagation of particle hole(ph) as well as delta- 
hole(A/i) excitations |43j. These considerations lead to modified hadronic tensor components (Jrpa) involving the 
bilinear terms, for which expressions are given in Refs. [13], [HI- 

III. INELASTIC RESONANCE PRODUCTION OF PIONS 

A. INCOHERENT PION PRODUCTION 

The basic reaction for the inelastic one pion production in nuclei, for a neutrino interacting with a nucleon N inside 
a nuclear target is given by 



v^k) + N(p) -» p-(fc') + N'[p') + 7T + (k„). 



(5) 



In the case of incoherent pion production in A dominance model, the weak hadronic currents interacting with the 
nucleons in the nuclear medium excite a A resonance which decays into pions and nucleons. The pions interact with 
the nucleus inside the nuclear medium before coming out. The final state interaction of pions leading to elastic, charge 
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FIG. 5: Ratio of the cross sections for induced charged current one pion production process and charged current quasielastic 
process in polystyrene. Experimental points have been taken from Ref.[l[ for K2K results in Scibar and from Refs.[48(|. (4i| for 
ANL results. 
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FIG. 7: Muon momentum distribution for induced charged current quasielastic process in 12 C averaged over the K2K 
spectrum. 
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FIG. 8: Q 2 distribution for induced incoherent charged current one pion production process in 12 C averaged over the K2K 
spectrum. 




FIG. 9: Muon momentum distribution for induced incoherent charged current one pion production process in C averaged 
over the K2K spectrum. 



exchange scattering and the absorption of pions lead to reduction of pion yield. The nuclear medium effects on A 
properties lead to modification in its mass and width which have been discussed earlier by Oset et al. pij and applied 
to explain the pion and electron induced pion production processes from nuclei. 

In the local density approximation the expression for the total cross section for the charged current one pion 
production (kinematics is shown in Fig.l) is written as 

i r~v ^ , f Q -^ r k '-^ r\^ 



(4..--.A ..... {Pp(r) + 9 Pn{r)W J Q ^ dQ . . ' /;/ ./ , dcosd " ./„ — Mm:. 

EEi^i 2 (6) 



E' + Err ( 1 — M-COS 



P 



where the proton density p p (r) = \pW) and the neutron density p n (r) = A ^ p(f) with p[r) as the nuclear density 
taken as 3-parameter Fermi Density [451 ]. 

The transition matrix element M. fi is given by 



V% m 

where L a is the leptonic current defined by Eq.([3]), and 0^ a is the N — A transition operator taken from Lalakulich 



M fi = ^S(p')^P ff ,0 Aa L ttU (p) (7) 

v 2 Tn-n 
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FIG. 10: Pion angular distribution for f M induced incoherent charged current one pion production process in C averaged over 
the K2K spectrum. 




FIG. 11: Pion Momentum distribution for induced incoherent charged current one pion production process in C averaged 
over the K2K spectrum. 



et al. 46|. 9\y is the weak mixing angle. V aX is the A propagator in momentum space and is given as : 

v ° = P 2 - M 2 A + zM A r (8) 

where P aX is the spin-3/2 projection operator and the delta decay width T is taken to be an energy dependent P-wave 
decay width Q: 

rW = -i(^) 2 §| qcm | 3 (9) 

67T \ TOtt / W 

1 1cm | is the pion momentum in the rest frame of the resonance and W is the center of mass energy. Inside the nuclear 
medium, the mass and width of delta are modified which in the present calculation are taken by considering the 
following effects. 

In nuclear medium As decay mainly through the A — > Nn channel. The final nucleons have to be above the Fermi 
momentum pp of the nucleon in the nucleus which leads to a modification in the decay width of delta. The modified 
delta decay width f has been taken from the works of Oset et al. [IJ which is given by 



f = Tx F{ PF ,E Al k A ) 



(10) 
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where 



F( PF ,E A ,k A ) = (ID 



E F = ^JM 2 + p' 2 F , k A is the A momentum and E A = yJW + k' A . 

Furthermore, in the nuclear medium there are additional decay channels open due to two and three body absorption 
processes like AiV — * NN and ANN — > NNN through which A disappears in nuclear medium without producing a 
pion, while a two body A absorption process like AN — > irNN gives rise to some more pions. These nuclear medium 
effects on the A propagation are included by modifying the mass and the decay width of A in nuclear medium as 

r^f-2im£ A and Ma -> Ma = Ma + ReE A . (12) 
The expressions for the real and imaginary part of the A self energy are taken from Oset et al. [44| : 

i?e£ A = iO—MeV and 
Po 

P \ , ^ P V , /i I P 



+ + C„ [f J (13) 

In the above equation Cq accounts for the AN — > ttNN process, Cai for the two-body absorption process AiV — > NN 
and Ca3 for the three-body absorption process ANN — > NNN. The coefficients Cq, Ca2, Ca3 an d a, (3 and 7 are 
taken from Ref. [i3| . 

The pions produced in this process are scattered and absorbed in the nuclear medium. This is treated in a Monte 
Carlo simulation which has been taken from Ref. [39(. 

B. COHERENT PION PRODUCTION 

f M induced coherent one pion production on 12 C target is given by +g 2 C — > /i~ +g 2 C + 7r + for which the basic 
reaction is given by Eq.([5]). The total cross section (kinematics is shown in Fig.l) is given by 

- — r Qma *dQ* f^dk' C dco^ d^&A ^\M fl \ 2 (14) 



(27T) 



where the matrix element for the reaction is given by 

L^ t is the leptonic current given by Eq.(3) and J M is the hadronic current given by 



M f i = -jpCosOcLUJuFtf- K) (15) 



j = V3^F(F 2 ) V u s (p')K a V° x 0^u r {v) (16) 



where F(P 2 ) (P=p+q) is ANn form factor given by[47j: 



- A 4 
F(P) ~ A 4 + (P 2 - M 2 ) 2 



with A = lGeV. 

J-(q — fcjr) is the nuclear form factor, given by 



P P {r) + \pn{r) 



When pion absorption effects are taken into account using the Eikonal approximation then the nuclear form factor 
!F(q — k n ) is modified to !F(q — k^), which is taken to be [37l ]: 

/>00 />00 

- M = 2tt / dz p(6,z) J (fe*6) e i(M-fcl)^ e -i/(M) (19) 



I) 
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where 

r°° 1 

f(b,z) = / -^Mp&z'^dz' (20) 
Jz 2|fc, r | 

A4 and fc£. are the longitudinal and transverse component of the pion momentum and the pion self-energy II is given 
by 

n(KM')) = | f^Y §Jki 2 p(M — gr (2i) 

with W as the center of mass energy in the A rest frame. 



IV. RESULTS AND DISCUSSIONS 



A. TOTAL CROSS SECTIONS 



For the quasielastic reaction, the numerical results are obtained from Eq.([2]) using vector form factors Fi(Q 2 ) and 
-^MQ 2 ) given by Bradford et al. [42j with vector dipole mass My=0.84GeV and a dipole form for Fa(Q 2 ) with axial 
dipole mass Myi=l.lGeV, and are shown in Fig. 2. Our results for the total scattering cross section a without the 
RPA effects (dotted line) in the case of charged current quasielastic scattering is in fair agreement with the calculation 
done in Fermi Gas Model Q > Hi > @ hi the energy region of 0.4-3 GeV. When RPA correlations are included in our 
model, we see that they give rise to a reduction in the total cross section a(E) (solid line) which is 12-15% in the 
energy region of 0.4-3 GeV. 

In the case of charged current induced incoherent and coherent pion productions the results for the total cross 
sections are obtained from Eq.© and Eq. (fT4]) respectively using the N-A transition form factors given by Lalakulich 
et al.4f| with Ma=1.1 GeV. The numerical results are shown in Fig. 3 and Fig. 4 for the incoherent and coherent 
pion productions. We see that the effect of nuclear medium in the pion production process as well as the final state 
interaction of pions both reduce the cross section. In the case of incoherent production of pions, the reduction due 
to nuclear medium effects in the production process is larger than due to final state interaction while in the case 
of coherent pion production, the reduction due to final state interaction is quite large as compared to the reduction 
due to the nuclear medium effects. However, the contribution of coherent process to the total pion production is 
small(5-6%) in the energy region of OAGeV < E < 3GeV, therefore, the reduction in the total cross section is mainly 
given by the reduction in the cross section of the incoherent process. 

In Fig. 5, we show the results for the ratio of total one pion production and the quasielastic production cross sections 

for induced reaction in polystyrene(Cs-f/8)j i- e - R(E) = a acQE( E { (solid line) as a function of neutrino energy. A 
comparison with the experimental results of K2K [l[ and ANL [Hj], [4{| for R(E) are also shown in the figure. In this 
figure, we also show the ratio R(E) (dashed line) when no nuclear medium effect is taken into account either in the 
quasielastic process or in the inelastic process. When the cross sections are calculated without any nuclear medium 
modifications except the Fermi motion and Pauli principle for the quasielastic process and without any modification 
of A properties in the nuclear medium in case of inelastic process, the results for R(E) are shown by dashed-double 
dotted line. This should be compared with the Monte Carlo predictions used in K2K analysis[l|. We see from Figs. 
2-5 that the nuclear medium modifications play an important role in quasielastic and inelastic reactions and tend to 
reduce the cross section ratio R(E). 

We have also studied the effect of the uncertainties in the ratio R(E) due to the use of various vector form factors 
given by Budd et al. [50| , Bosted [5l[ and dipole form factors (53 in the quasielastic case and the N-A transition form 
factors given by Schreiner & von Hippel[53[ and Paschos et al. [32j| in the inelastic case. The effect of varying the axial 
dipole mass (Ma) in the parameterization of axial form factor in quasielastic as well as in the inelastic processes in 
the region of 1.05 < Ma < 1.21 GeV on the ratio R(E) has also been studied. These effects lead to a small change of 
about 4-6% in the ratio R(E). 



B. ANGULAR AND MOMENTUM DISTRIBUTIONS OF LEPTONS AND PIONS 



In Figs. 6 & 7, we show the effect of nuclear medium on Q 2 distribution and the momentum distribution of the 
leptons averaged over the neutrino flux at K2K in the case of neutrino induced charged current quasielastic scattering 
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on 12 C target. A reduction in Q 2 distribution in the peak region is found to be about 30% when RPA effects are 
taken into account while this reduction is about 12% in the momentum distribution of leptons in the peak region. 
In the case of incoherent one pion production process, the Q 2 distribution as well as the momentum distribution 
of leptons are shown in Figs. 8 & 9. In these reactions charge pions are also produced and we show in Figs. 10 & 
11, the angular distribution as well as momentum distribution of pions and the effect of nuclear medium on these 
distributions. We find that the reduction in Q 2 distribution in the peak region is around 30% when nuclear medium 
effects arc taken into account which further reduces by about 14% when pion absorption effects are incorporated. For 
the pion momentum distribution the reduction due to nuclear medium effects is around 40% which further reduces 
by about 15% when pion absorption effects are taken into account. We have not shown these results for the coherent 
process as the contribution to the total pion production events have been found to be quite small. However, these 
results are interesting in their own right and will be reported elsewhere. 

To conclude, we have studied in this paper the effect of nuclear medium on the inclusive quasielastic and one pion 
charged current production in nuclei and compared our results with the recent experimental results from the K2K 
experiment [l| . The nuclear medium effects are important in this energy region and reduce the ratio by 30% around 
the energy region of lGeV which becomes 20-25% in the neutrino energy region of 2-3GeV. The results for nuclear 
medium effects on the Q 2 and the momentum distribution of lepton in the quasielastic as well as inelastic reaction 
and the results for the nuclear medium effects on the angular and momentum distribution of pions have also been 
presented. 
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